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ABSTRACT 

This paper describes the processing applied to the Planck High Frequency Instrument (HFI) cleaned, time-ordered information to produce pho¬ 
tometrically calibrated maps in temperature and (for the first time) in polarization. The data from the entire 2.5 year HFI mission include almost 
five independent full-sky surveys. HFI observes the sky over a broad range of frequencies, from 100 to 857 GHz. To obtain the best accuracy on 
the calibration over such a large range, two different photometric calibration schemes have been used. The 545 and 857 GHz data are calibrated 
using models of planetary atmospheric emission. The lower frequencies (from 100 to 353 GHz) are calibrated using the time-variable cosmological 
microwave background dipole, which we call the ’’orbital dipole”. This source of calibration only depends on the satellite velocity with respect to 
the solar system. Using a CMB temperature of Tcmb - 2 .7255 + 0.0006 K, it permits an independent measurement of the amplitude of the CMB 
solar dipole (3364.3 + 1.5 pK) which is approximatively Icr higher than the WMAP measurement with a direction that is consistent between both 
experiments. 

We describe the pipeline used to produce the maps of intensity and linear polarization from the HFI timelines, and the scheme used to set the zero 
level of the maps a posteriori. We also summarize the noise characteristics of the HFI maps in the 2015 Planck data release and present some null 
tests to assess their quality. Finally, we discuss the major systematic effects and in particular the leakage induced by flux mismatch between the 
detectors that leads to spurious polarization signal. 

Key words, cosmology: observations - cosmic background radiation - surveys - methods: data analysis 
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1. Introduction 

This paper, one of a set associated with the 2015 Planck data 
release, is the second of two describing the processing of the 
data from the High Frequency Instrument (HFI). The HFI is 
one of the two instruments on board Planck^, the European 
Space Agency’s satellite dedicated to precision measurements of 
the cosmic microwave background (CMB). The HFI uses cold 
optics (at 4K, 1.6 K, and 100 mK), filters, and 52 bolometers 
cooled to 100 mK. Coupled to the Planck telescope, it enables 
us to map the continuum emission of the sky in intensity and po¬ 
larization at frequencies of 100, 143, 217, and 353 GHz, and in 
intensity at 545 and 857 GHz. Paper 1 (Planck Collaboration VII 
2015) describes the processing of the data at the time-ordered 
level and the measurement of the beam. Paper 2 (this paper) de¬ 
scribes the HFI photometric calibration and mapmaking. 

Paper 1 describes how the TOI (time-ordered information) 
of each of the 52 bolometers is processed and fiagged. Sampled 
at 5.544 ms, the TOI is first corrected for the ADC nonlinearity, 
then it is demodulated and converted to the absorbed power with 
a simple nonlinear bolometric correction. Glitches are fiagged 
and glitch tails are removed from the TOI. Thermal fiuctuations 
are removed on the 1 min timescale. Sharp lines in the tempo¬ 
ral power spectrum of the TOI from the infiuence of the 4-K 
cooler are removed. Finally, the bolometer time response is de¬ 
convolved and the TOI is low-pass filtered. At this point, the 
TOIs are cleaned but not yet calibrated. The measurement of the 
beam is done using a combination of observations of planets for 
the main beam and GRASP calculations for the sidelobes. The 
focal plane geometry, or the relative position of bolometers in 
the sky, is deduced from Mars observations. 

This paper describes how the prepared TOIs are used to make 
the calibrated maps for all Planck HFI bands. After a summary 
of the photometric definitions (Sect. 2), this paper gives a de¬ 
scription of the main steps of the mapmaking processing, fo¬ 
cusing on the changes made since the previous Planck releases 
(Planck Collaboration VI 2014; Planck HFI Core Team 2011). 
The major difference concerns the calibration (Sect. 3) which 
is now based on the orbital CMB dipole for the lower-frequency 
channels (100, 143, 217, and 353 GHz, also called ”CMB chan¬ 
nels”) while the 545 and 857 GHz channels are photometrically 
calibrated using the signal from Uranus and Neptune. Section 4 
explains the mapmaking upgrades, including the polarization 
treatment. Section 5 describes the maps, the solar dipole mea¬ 
surement, and the derivation of far sidelobes and zodiacal maps. 
Section 6 presents the noise characteristics and the null tests ob¬ 
tained by splitting the Planck HFI dataset into different groups 
based on ring period, time period, or detector sets. Consistency 
checks are performed in order to assess the fidelity of the maps. 
Finally, Sect. 7 is dedicated to the description of systematic ef¬ 
fects, in particular in polarization. The major systematic resid¬ 
uals in Planck HFI data are due to the leakage from tempera¬ 
ture to polarization induced by flux mismatch between associ¬ 
ated bolometers. This is the result of either bandpass mismatch, 
or zero level uncertainty, or calibration uncertainty. We present a 
first attempt to correct the maps for the intensity-to-polarization 
leakage. At lower frequencies, even after correction, the residu- 


^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci¬ 
entific consortia funded by ESA member states and led by Principal 
Investigators from France and Italy, telescope reflectors provided 
through a collaboration between ESA and a scientific consortium led 
and funded by Denmark, and additional contributions from NASA 
(USA). 


als are still higher than the noise level and thus the maps cannot 
yet be used for cosmological analysis and are not released. 

2. Photometric Equations 

The power absorbed by a given detector at time t can be written 
as a sum of three terms corresponding to the first three Stokes 
parameters (7^, Qp, Up) at the sampled pixel p of the beam- 
convolved sky: 

Pf — G -l- p ^Qp cos F tr) U p sin + ^)|] F , (1) 

where G encodes a photometric calibration factor, p is the de¬ 
tector polarization efficiency, if/t is the roll angle of the satel¬ 
lite, a is the detector polarization angle, and nt represents all 
the noise contributions to the absorbed power (photon noise, 
phonon noise, glitch residuals, etc.). The polarization efficiency 
is derived from the cross-polarization coefficient p through p = 
(1 - 77 )/(1 -r 77 ). It allows us to describe spider-web bolometers 
(SWB, p 0) as well as polarization-sensitive bolometers (PSB, 
p ^ 1). According to the bolometer model and given the stability 
of the HFI operational conditions during the mission, the gain G 
is expected to be constant over the whole mission (see section 1 
of Paper 1), once the bolometer nonlinearities have been cor¬ 
rected. 

For an axisymmetric beam response, the “smearing” and 
“pointing” operations commute, and one can solve directly for 
the pixelized beam-convolved map: 

Pt = G(l,pcos2(i//t-\-a),psm2(i//t-\-a))'^ •(^Ip,Qp,Up^-\-nt 

= GxAtpTp -\-nt, ( 2 ) 

where the direction of observation at the time t, (6t, (pt), falls into 
pixel p and we define the map-pointing matrix A and the sky 
signal T = (7, Q, U). 

3. Calibration 

The bolometer signal measured through current-biasing is pro¬ 
portional to the small variation in the incoming power from the 
sky. To express the measurement in sky temperature units, one 
has to determine a gain per detector based on a known source in 
the sky. For the HFI low frequency channels (100 to 353 GHz), 
we use the CMB orbital dipole as a primary calibrator. This sig¬ 
nal fills the entire beam and is almost insensitive to the beam 
profile and only marginally affected by pointing errors, while its 
signal-to-noise is high enough thanks to the full-sky coverage. 
Moreover, it is a stronger signal than CMB anisotropies (by a 
factor of around 10), but not bright enough to cause nonlineari¬ 
ties in the detectors, and has the same electromagnetic spectrum 
as the anisotropies. At higher frequencies (545 and 857 GHz), 
calibration is performed on planets. 

3.1. CMB dipole conventions 

The CMB dipole is induced by the Doppler effect of the relative 
motion of the satellite with respect to the CMB frame, 

T'DopplerCf,«) = .. , (3) 

7^(1-A •«) 

where Pt = Vtic and yt = (1 - Vt is the satellite velocity 

at time t and u is the unit vector along the line of sight. 
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The solar system motion with respect to the CMB frame, giv¬ 
ing rise to what is referred to as the “solar dipole”, is the domi¬ 
nant component of the satellite velocity. A residual contribution 
(called the“orbital dipole”) is induced by the yearly motion of 
the satellite with respect to the solar system barycentre. The so¬ 
lar dipole can be considered as sky-stationary during the obser¬ 
vations and is thus projected onto the sky as an ^ = 1 compo¬ 
nent with amplitude previously measured by COBE and WMAP, 
3355 + 8 yuK (Hinshaw et al. 2009). Relativistic corrections to 
the solar dipole produce second order anisotropies at multipoles 
^ > 1, with amplitudes proportional to , and more importantly 
couple the two dipole components, as will be discussed below. 
Though the orbital dipole velocity is typically an order of mag¬ 
nitude lower than the solar dipole, it is time dependent, and its 
time-variability is precisely determined by the satellite velocity 
which is known at the level of 10“"^ km s“^ Finally, to calibrate 
in temperature, we only rely on an external measurement of the 
CMB absolute temperature. We use Tqmb = 2.7255 K (Fixsen 
2009). 

The expansion of Eq. (3) in yS gives 

^ CMB ^ 

If we decompose the velocity into a solar boost Pi and an orbital 
boost P2 then 

A = Al + A2 -ySi -ySz + 2081 • u){fi2 ■ u) + 0 (J 3 ^) ( 5 ) 

where the first term corresponds to the solar dipole, the second 
term is the orbital dipole, and the third and fourth terms show 
the coupling between both due to relativistic corrections. 

3.2. Absolute calibration on orbital dipole 

The calibration algorithm takes advantage of the orbital dipole 
not being fixed on the sky, unlike the solar dipole (during the 
length of the mission). In practice, relativistic corrections and 
second order in the development of the conversion from ly 
to Tcmb couple solar and orbital dipoles creating an additional 
non-stationary signal which also depends on the frequency. We 
use as calibration reference signal, the total CMB dipole com¬ 
puted using Eq. (3), assuming that it has the same scaling with 
frequency than the higher multipole CMB anisotropies. This 
approximation will leave in the HFI frequency maps the fre¬ 
quency dependent fraction of the kinematic quadrupole aris¬ 
ing from the 2d order term in the p expansion described in 
(Kamionkowski & Knox 2003) and (Quartin & Notari 2015). 
The amplitude of the kinematic quadrupole is expected to be 
lower than 0.5% on the CMB-calibrated PlanckAAFl channels. 
However, it is not directly correlated to the orbital dipole on 
which we calibrate the data. Thus we expect the systematic error 
induced on the gain estimation to be much smaller. 

The orbital dipole signal is modulated on a one-year period. 
To take into account the time variation of this signal, we need to 
add a term in Eq. (2): 

P = G x. (AT -i- ^orb) F ri , (6) 

where ^orb is the time-dependent orbital dipole signal, while in 
this formula the solar dipole is part of the sky signal T. Note that 
we can arbitrarily set all or part of the solar dipole signal either 
in the calibration template or in the sky without changing the 
resulting gain estimation. 


Since ^orb is known, we can solve Eq. (6) for each bolometer 
independently, rewriting the system as 

P =■ AT -b C^orb F ri , (7) 

where the unknowns are the sky-signal in absorbed power units 
T and the gain G. The calibration problem is thus linear and can 
be solved directly. The maximum-likelihood estimate of the gain 
G is obtained by combining all available samples using the noise 
covariance matrix N = and marginalizing over T, solving 
the equation 

(<lN-iZVb)xG = <lN-iZrf, (8) 

where d is the vector of input data Pf and we define 

Z= I-A(aTN-‘a)“‘A^N- b (9) 

In practice, the noise is treated by assuming that the de- 
striping reduces the matrix N to a diagonal one once the data d 
have been corrected for a constant offset (see Sect. 4.2). For the 
HFI, the mapmaking problem is degenerate for the reconstruc¬ 
tion of polarization if we solve for a single detector at a time. 
Neglecting polarization in Eq. (8) for polarization-sensitive de¬ 
tectors biases the calibration solution. Moreover, we need a very 
accurate relative calibration between detectors that are combined 
to reconstruct polarization in order to minimize leakage from in¬ 
tensity to polarization (see Sect. 7.3). For this reason, for the 
Planck 2015 release we have extended the algorithm described 
in Tristram et al. (2011) to perform a multi-detector gain esti¬ 
mation for all bolometers at a given frequency together with the 
offsets (see Sect 4.2). 

3.3. Long time-constant residuals 

The observation strategy of Planck results in the path across 
a particular part of the sky being almost reversed 6 months 
later. As described in Paper 1 , we take advantage of this to de¬ 
rive the time transfer function below one second. Nevertheless, 
longer time responses (larger than the second), even with low 
amplitudes, may bias the calibration estimation by distorting the 
dipole signal and causing some leakage of the solar into the or¬ 
bital dipole. To take into account the systematic residuals af¬ 
ter time-constant deconvolution, we built a simplified model de¬ 
scribing a pure single-mode sinusoidal dipole signal (including 
solar and orbital dipole) convolved with an exponential decay in 
the time domain. In the frequency domain, this reads: 

Tdip = r(tdip * Be *''') = ’ (10) 

where B and r are the amplitude and the time constant, and Vspin 
is the spin frequency. Tdip is a complex coefficient, the real part 
of which corresponds to the relative change in the gain G, while 
the imaginary part corresponds to the amplitude of the dipole 
mode shifted by 90°. 

In practice, prior to the absolute calibration, we solve for 
the amplitude of a shifted-dipole template using single detector 
calibration: 

P = Af + G/dip+O^fJ + n, (11) 

where C gives the amplitude of the dipole mode shifted by 90°. 
With this toy-model, we cannot reconstruct the amplitude and 
the time constant because we only fit for one coefficient C. But 
we trace the systematic effect on the dipoles due to the very long 
time response. Once the coefficient C has been determined, we 
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correct the data for to account for this additional shifted 
mode coming from the residual time constant. 

The effect on the gain G depends on the unknown value of 
the time constant r. Figure 1 shows the impact on the gains for 
each bolometer when including the shifted-dipole correction. At 
higher frequencies (353 GHz and above), the signal is no longer 
dominated by a dipole and cannot be approximated by the model 
described above. For the 2015 release, we do not include any cor¬ 
rection for time constant residuals in the mapmaking for those 
channels. 


0.3 
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Fig. 1. Effect on the gains when including the model for long 
time constant residuals as described in the text. 
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3.4. Submillimetre calibration 

Mars, Jupiter, Saturn, Uranus, and Neptune are all observed sev¬ 
eral times during the full length of the mission. The submillime¬ 
tre channels of HFI (545 and 857 GHz) are calibrated on models 
of Uranus and Neptune. We do not use Jupiter and Saturn obser¬ 
vations for calibration, since both planets have strong absorption 
features at those frequencies, which complicate comparison with 
broadband measurements. The flux from Jupiter also leads to de¬ 
tector saturation at the highest HFI frequencies. Similarly, we 
choose not to use Mars as a calibrator, because strong seasonal 
variations complicate the modelling. 

Various methods are used to derive planet flux densi¬ 
ties, including aperture photometry and point-spread function 
(PSF) fits. Planet measurements with Planck are studied fully 
in Planck Collaboration XXXIV (2015). We focus here on cali¬ 
bration using aperture photometry in the time-ordered data from 
the submillimetre channels. The simulation pipeline used for the 
main beam reconstruction computes the reconstruction bias and 
the error budget (the mean and the variance) for each planet ob¬ 
servation. The comparison of flux measurements of Neptune and 
Uranus with up-to-date models provides the calibration factors 
at 545 and 857 GHz. 

3.4.1. Planet models 

We use the ESA 2 model for Uranus and the ESA 3 model for 
Neptune (Moreno 2010). Both models quote absolute uncertain¬ 
ties of 5 %. Planet model spectra are produced from their mod¬ 


elled brightness temperatures using the planet solid angles at 
the time of observation and integrated under the individual 545 
and 857 GHz bolometer bandpasses. Flux densities are colour- 
corrected to a reference spectrum defined by a constant v/y law, 
so as to be directly comparable to HFI flux density measure¬ 
ments. 

Planck Collaboration XXXIV (2015) give a detailed account 
of the ratio between the expected planet fluxes and the mea¬ 
sured ones at all HFI frequencies and for the five observed plan¬ 
ets (Mars, Jupiter, Saturn, Uranus, and Neptune). They all fall 
within the 5 % model uncertainty range. This is a validation of 
the models at 100-353 GHz. Hence the models can be used with 
some confidence to calibrate the 545 and 857 GHz channels. 

3.4.2. Aperture photometry in the timelines 

We select all samples in a 2° x 2° box around the planet positions 
and build time-ordered vector objects (hereafter “timelines”) for 
each bolometer and each planet scan. We use the first four scans 
of Neptune and Uranus (season 1 to 4). We build correspond¬ 
ing background timelines using all the samples in a 2° x 2° box 
around the planet position when the planet is not there. The re¬ 
sulting background has a much higher spatial density than the 
planet timelines. We use this to build a background mini-map 
with 2' X 2' pixels that can then be interpolated at each sample 
position in the planet timelines in order to remove a local back¬ 
ground estimate (see appendix B in Paper 1 for details). 

The aperture photometry measurement procedure applied to 
our planet timelines is an extension of the usual aperture pho¬ 
tometry approach to irregularly gridded data. Flux is integrated 
in an aperture of radius 3 times the effective FWHM. Typically, 
aperture photometry is applied to maps of fixed-size pixels, 
which means integrating the flux in the aperture is equivalent to 
summing the pixel values. In our case, we have to take into ac¬ 
count the inhomogenous spatial distribution of the samples. To 
do this we assume that the beam products perfectly describe the 
spatial light profile of Neptune and Uranus, and compute a spa¬ 
tial sampling correction factor as the ratio of the integrated flux 
in a highly spatially oversampled beam and the integrated flux in 
a beam sampled on the planet timelines. This sampling correc¬ 
tion has to be computed for each bolometer and planet crossing, 
because the spatial sampling varies between planet observations. 
We estimate the statistical uncertainty of the measurements as 
the standard deviation of the samples in an annulus of radius 3 
to 5 times the effective FWHM of the beam. 

We And large variations between the individual planet mea¬ 
surements in each detector and at each season of the full mis¬ 
sion survey (Fig. 2) which we attribute to underestimation of 
the measurement uncertainty. Indeed, the signal from Neptune 
and Uranus is not expected to vary in time apart from the dif¬ 
ferences in solid angle which are very small and already taken 
into account. While accurately corrected by the timeline aper¬ 
ture photometry algorithm presented here, the limited available 
spatial sampling of the planet signal at these frequencies could 
explain part of the variations. We therefore decided to include 
the “seasonal” rms of the measurements (we have 4 observations 
per bolometer per planet) in the measurement uncertainty. 

The averaged calibration factors for each detector for each 
planet are in very good agreement. The final calibration factors 
are the means of both planet estimates. We compare them to the 
2013 Planck release in Fig. 3: the calibration factors changed 
by 1.9 and 4.1% at 545 and 857 GHz, respectively, which is 
within the planet modelling uncertainty. Combined with other 
pipeline changes (such as the ADC corrections), the 2015 fre- 
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Fig. 2. Dispersion of the planet-derived calibration factor per 
season around the planet calibration estimates for Uranus (top), 
and Neptune (bottom). 



Fig. 3. Comparison of the 2015 Planck release calibration to the 
one from the 2013 release. We show the relative difference in 
percent per bolometer for Uranus (green), Neptune (blue), and 
both calibrators combined (black). 


quency maps have decreased in brightness by 1.8 and 3.3 % com¬ 
pared to 2013. 

Calibration uncertainties are given in Table 1 . In order to pro¬ 
duce the frequency maps, detectors are weighted by their inverse 
noise variance. We use the same weights to compute the cor¬ 
responding calibration errors. We estimate combined statistical 
errors of 1.1 % and 1.4% at 545 and 857 GHz respectively, to 


Table 1. Uncertainties in the planet-derived calibration factor for 
each bolometer. The systematic uncertainty is the absolute un¬ 
certainty of the planet model. 


Bolometer Uncertainty Uncertainty 
Stat. Syst. 

545-1 . 1.0% 5.0% 

545-2 . 1.8% 5.0% 

545-4 . 2.3% 5.0% 

857-1 . 2.6% 5.0% 

857-2 . 2.9% 5.0% 

857-3 . 2.8% 5.0% 

857-4 . 2.0% 5.0% 


which one should add linearly (as it should be done for system- 
atics), the 5 % uncertainty arising from the planet models. Errors 
on absolute calibration are therefore 6.1 and 6.4% at 545 and 
857 GHz, respectively. Since the reported relative uncertainty of 
the models is of the order of 2 %, we find the relative calibration 
between the two HFI highest frequency channels to be better 
than 3 %. 

4. Mapmaking 

4.1. Summary 

Each data sample is calibrated in Kcmb for the 100,143, 217, and 
353 GHz channels, or MJy sr“^ (assuming a constant vly law) for 
the 545 and 857 GHz channels, using the calibration scheme de¬ 
scribed above. Unlike for the 2013 release, the bolometer gains 
are assumed to be constant throughout the mission. The Planck 
total dipole (solar and orbital) is computed and subtracted from 
the data. 

As in the Planck 2013 release, we average the measurements 
for each detector in each pixel visited during a stable point¬ 
ing period (hereafter called a “ring”) while keeping track of 
the bolometer orientation on the sky. The subsequent calibra¬ 
tion and mapmaking operations use this intermediate product as 
input. The calibrated TOIs are only modified by a single offset 
value per ring and per detector. The offsets are determined with 
a destriping method described in Tristram et al. (2011). Here, 
the size of the pixels where the consistency of different rings is 
tested is 8' (A^side = 512). Maps in intensity and polarization are 
used to assess the consistency of the destriper solution. The off¬ 
sets are simultaneously determined for all bolometers at a given 
frequency using the full mission data. For all the maps produced 
using a given bolometer, the same offset per ring is used (except 
in the case of half-rings, see below). 

The products of the HFI mapmaking pipelines are maps of 
(/, 2 , U) together with their covariances (//, IQ, lU, QQ, QU, 
UU), pixelized according to the HEALPix scheme (Gorski et al. 
2005). The map resolution is Aside = 2048, and the pixel size 
is 1.'7. The mapmaking method is a simple projection of each 
unfiagged sample to the nearest grid pixel. In the case of polar¬ 
ization data with several detectors solved simultaneously, the po¬ 
larization equation is inverted on a per-pixel basis (see Sect. 4.2). 
For each sky map, a hit count map is computed (number of sam¬ 
ples per pixel; one sample has a duration of 5.544 ms). 

We provide maps from which the zodiacal light compo¬ 
nent, which varies in time, has been removed, based on tem¬ 
plates fitted on the survey-difference maps (see Sect. 5.4). These 
templates are systematically subtracted from the data of each 
bolometer prior to the mapmaking. 
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Unlike in the 2013 release, the far sidelobes (FSL) are not 
removed from the maps. At most, this leaves residuals of or¬ 
der 0.5-1.5 yuK in the 100-353 GHz maps, with uncertainties 
on the residuals of roughly 100 % (Planck Collaboration VII 
2015). At higher frequencies. Galactic pick-up from the FSL 
produces significant residuals in about half of the detectors of 
the order of 0.03 MJy sr“^ at 545 and 0.3MJysr“^ at 857 GHz 
(see Planck Collaboration X 2015). The total solid angle in the 
spillover is a better known quantity, and we describe in Sect. 5.3 
the effect of the FSL on the calibration. 

4.2. Mapmaking method 

In the same way as in Planck Collaboration VIII (2014), we use 
a destriping algorithm to deal with the HFI low-frequency noise. 
In this approach, the noise in a ring r is represented by an offset, 
denoted Or, and a white noise part n, which is uncorrelated with 
the low-frequency noise. For a given bolometer, we can write 
Eq. (6) as: 

Pt — G (Atp Tp ^orb) ^tr Or ■> (1^) 

where L^^ is the ring-pointing matrix (which associates the data 
sample t with the ring number r). The unknowns are the gain G, 
the offsets for each ring Or, and the sky signal for each pixel Tp = 
(Ip, Qp, Up). As there is a degeneracy between the average of the 
offsets and the zero level of the maps, we impose the constraint 
(Or) = 0. The absolute zero level of the maps is determined as 
described in Sect. 4.4. 

For the production of the maps for the 2015 HFI data release, 
we first build the rings for all detectors. We apply the following 
frequency-dependent processing to these compressed data sets. 
For CMB frequencies (100, 143, and 217 GHz) channels: 

1. we estimate a first approximation of the orbital dipole gain 
together with the offsets and the amplitude of the long¬ 
time-constant residuals for each bolometer independently, 
neglecting the polarization signal; 

2. we then derive the gains and offsets for a fixed amplitude of 
the long-time-constant residuals using the multi-bolometer 
algorithm. 

At 353 GHz, the long-time-constant residuals are more dif¬ 
ficult to constrain. They are driven more by Galactic emission 
drifts than by the dipole, which dominates at lower frequencies, 
so that the model described in Sect. 3.3 is not relevant. Hence 
for this frequency we use a simpler pipeline without a long-time- 
constant residuals template: 

1. we estimate the orbital dipole gain together with the offsets 
for each bolometer independently, neglecting the polariza¬ 
tion signal; 

2. we then estimate the final offsets using a destriping proce¬ 
dure for all bolometers at this frequency. 

For the two highest frequencies, at 545 and 857 GHz, the 
pipeline is considerably different, because we use the planets 
(Uranus and Neptune) as calibration sources: 

1. we estimate a first approximation of the offsets for each 
bolometer independently; 

2. we derive the gains from the planet flux comparison; 

3. we then estimate the final offsets with a destriping procedure 
for all bolometers at a given frequency. 


Finally, for each data set, using the pre-computed gains and 
offsets, we project the ring data onto maps. For polarization, we 
invert the 3x3 system derived from Eq. (12) for each pixel in¬ 
dependently with a criterion that the condition number be lower 
than 10^: 

r= (A^N“^a)”‘(13) 

where d are the calibrated, cleaned ring data (offsets and or¬ 
bital dipole removed) d = (P -V o) IG - Grb- Note that we use 
HEALPix (not lAU) conventions for the sign and normalization 
of the Q and U Stokes parameters. 

For destriping, we use the same tool as before, polkapix, 
which was thoroughly validated in Tristram et al. (2011). Maps 
are built by simple co-addition in each pixel of the destriped, 
calibrated, and time-varying component-subtracted signal. We 
subtract the CMB dipole as measured by Planck (see Sect. 5.1). 

We introduced the following modifications with respect 
to Planck Collaboration VIII (2014): 

- we include polarization in the destriping for the channels that 
include PSBs; 

- we enlarged the masked fraction of the sky, from 10 to 15 %, 
based on Galactic emission, to avoid signal gradients leaking 
into offsets; 

- to improve the offset accuracy, we compute one set of off¬ 
sets combining all detectors, using full-mission data, and use 
them for all derived maps involving these detectors. 

This last change induces a small noise correlation between 
detector-set maps (see Sect. 6.5). In 2013 we computed inde¬ 
pendent offsets for each detector or detector set, including the 
full mission. 

4.3. Map products 

The principal HFI final product consists of six maps that cover 
the six frequencies (100-857 GHz) for the full mission in in¬ 
tensity at high resolution (Aside = 2048). However, many more 
maps are needed to assess the noise and the consistency of 
the data. Figure 4 summarizes the various splits produced. The 
branches that are delivered in the 2015 release are described in 
Appendix A. 

Maps from different halves of each ring period (first and last) 
are computed independently of each other, including the offset 
per half-ring. Thus, half-ring half-difference maps can give a 
quick account of the noise level in the maps. 

For each frequency, we also produce temperature and polar¬ 
ization maps using detector sets (each set including four polar¬ 
ization sensitive bolometers). In addition, we produce a temper¬ 
ature map for each spider-web bolometer (SWB). 

Planck's scanning strategy samples almost all the sky pixels 
every six months, with alternating scan directions in successive 
6-month periods. The full cold HFI mission encompasses five 
surveys, each covering a large fraction of the sky. Surveys 1- 
2 and 3-4 are paired to produce Year 1 and Year 2 maps (see 
Table A.2). Maps are produced for the full-mission dataset to¬ 
gether with the survey, year, and half-mission maps. With each 
map is associated a hit-count map and variance maps (II; and 
QQ, UU, IQ, lU, and QU when polarization is reconstructed). 
Overall, a total of about 6500 sky maps have been produced. We 
have used this data set to evaluate the performance of the photo¬ 
metric calibration by examining difference maps (see Sect. 6.2). 
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Fig. 4. Map matrix. The HFI maps are released in different flavours. Not all combinations are released but any map will correspond 
to a choice of lines in each box. The Frequency box is related to the use of all detectors at a given frequency (Channel), or individual 
bolometers or sets of detectors as deflned in Table A. 1 . The Ring box is a way of splitting (or not) the data in equal halves at the 
ring level. The Duration box indicates the different ways of splitting data between surveys, years, full or nominal mission, first half 
mission, or last half mission. The Component box indicates the systematic corrections that can be applied at the map level. The 
recommended first choice map is highlighted in red. See Appendix A for details. 


4.4. Zero levels 

Planck-HFl cannot measure the absolute sky background. The 
mapmaking procedure does not change the mean value of 
the input TOI. We therefore adjust the monopole on the 
maps a posteriori, in a similar manner to the method used 
in Planck Collaboration VIII (2014) that relies on external 
datasets. To achieve this, we need to take into account two major 
components of the monopole: 

1. Galactic dust emission: we estimate the brightness in the HFI 
single-detector maps that corresponds to zero gas column- 
density (i.e., zero Galactic dust emission). To do so, we use 
the Hi column density, which is assumed to be a reliable 
tracer of the Galactic gas column-density in very diffuse ar¬ 
eas (see Planck Collaboration VIII 2014, Sect. 5.1). The off¬ 
sets derived are then subtracted from each detector’s data in 
the processing. 

2. Extragalactic emission: the cosmic infrared background 
(CIB) monopole is taken into account by adding the levels 
from Bethermin et al. (2012) to the maps (see Table 6). 

The sum of the two offsets is appropriate for total emission 
analysis. For Galactic studies, only the Galactic zero level has 
to be set which can be achieved by subtracting the CIB levels 
(Table 6) from the released maps. Unlike for the previous re¬ 
lease, in the 2015 maps the zero level correction (both CIB and 
Galactic) has been applied. 

Zodiacal light has not been accounted for in this procedure. 
The offsets that have to be removed at each frequency to set 
the Galactic zero level using zodiacal-light-corrected maps are 
smaller than those needed for the total maps. The correction to 
be applied to the released maps are given in Table 6. 

4.5. Polarization efficiency and orientation 

The calibration parameters for the PSBs were measured on the 
ground before launch. Rosset et al. (2010) have reported pre¬ 
flight measurements of the polarization efficiency of the HFI 
PSBs with an accuracy of 0.3 %. The absolute orientation of the 
focal plane has been measured at a level better than 0?3. The rel¬ 


ative orientation between PSBs is known with an accuracy better 
than 0?9. 

The SWBs are much less sensitive to polarization. 
Nonetheless, we take into account their polarization efficiency 
which is between 1 and 9 %, although their orientations have 
been less accurately determined (errors can be up to a few per¬ 
cent), as described in Rosset et al. (2010). 

5. HFI temperature and polarization maps 

5.1. Solar dipole measurement 

The €= I mode of CMB anisotropy is unique in that its am¬ 
plitude is dominated by a large component associated with our 
motion with respect to the CMB rest frame. In this section, we 
present the CMB solar dipole results based on Planck-HFl maps 
at the two lowest frequencies, 100 and 143 GHz. Low-frequency 
maps are dominated by CMB over a large fraction of the sky. 
Nevertheless, the inhomogeneous nature of the dust emission 
can bias CMB solar dipole estimates. 

We cleaned the Galactic emission from the HFI maps using 
a local correlation with the 857 GHz map. We model each HFI 
map /y as 

ly - C = qhsi + Aes (14) 

where C is the CMB anisotropy (here we use the SMICA 
map, Planck Collaboration IX 2015, from which we remove any 
residual dipole component) and l^si is the Planck 857 GHz map 
that is assumed to have a negligible contribution from the solar 
dipole.^ The term includes the dipole and any systematic 
effects from both I^si and ly. In bright regions of the sky, Dj-es 
also contains extra emission that is uncorrelated or only partially 
correlated with /gsv, for instance free-free emission or CO. 

In order to capture any spatial variations of the dust SED, 
we estimated q and Dj-es on an Aside = 64 grid. For each Aside = 
64 pixel, we performed a linear regression of the Aside = 2048 
pixels of ly vs / 857 , assuming a constant dust SED over a 55' 
area. We then fit for the dipole amplitude and direction in Dres 

^ The amplitude of the solar dipole at 857 GHz is 0.0076 MJy sr"^ 
At least 90 % of this is removed in the mapmaking process, leaving a 
residual that is well below the noise level and any systematic effects. 
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Table 2. CMB solar dipole measurements for the 100 and 143 GHz channels estimated for different sky coverage levels (37, 50, 
and 58%) corresponding to three thresholds in 857 GHz amplitude (2, 3, and 4 MJysr“^). Uncertainties include only statistical 
errors. Systematic errors are 0.8 yuK for the amplitude, and (0?024, 0?0034) in Galactic (longitude, latitude). 


Frequency Threshold d Ion lat 

[GHz] [MJysr-i] [yuK] n C] 


100 . 2 3364.81 + 0.06 263.921 + 0.002 48.2642 + 0.0008 

100 . 3 3364.76 + 0.05 263.922 + 0.002 48.2640 + 0.0006 

100 . 4 3364.99 + 0.04 263.928 + 0.002 48.2631 + 0.0006 

143 . 2 3364.05 + 0.03 263.908 + 0.001 48.2641 +0.0004 

143 . 3 3363.72 + 0.02 263.903 + 0.001 48.2653 + 0.0003 

143 . 4 3363.39 + 0.02 263.905 + 0.001 48.2668 + 0.0003 


using sky pixels where hsi < 2, 3, or 4 MJy sr"^ (corresponding 
to 37, 50, or 58 % of the sky respectively) to limit the effect of 
Galactic-emission residuals (CO, free-free emission, and small- 
scale dust SED variations). 

The results are given in Table 2. We measure a Solar System 
peculiar velocity of 370.06 + 0.09 km s“^ with respect to the 
CMB rest frame. We use the CMB temperature from Fixsen 
(2009) (2.7255 + 0.0006 K) to convert that measurement into a 
CMB dipole measurement. 

The error bars here only include statistical uncertainties, 
which are very low thanks to the Planck-W\ signal-to-noise 
ratio. We evaluate the additional systematic uncertainties from 
the variation of the results between independent bolometer 
maps. For the amplitude, the peak-to-peak variation between 
bolometers and combined maps is +0.8 yuK at 100 and 143 GHz. 
Variations with sky coverage are of the same order. Note that 
the uncertainty from the FIRAS temperature should be added to 
the budget (+0.74yuK). For the coordinates, we found variations 
of +0?013 in longitude and +0?0019 in latitude. These differ¬ 
ences are observed when comparing results at different frequen¬ 
cies, and are likely to result from uncertainties in the foreground 
subtraction. This is also consistent with the magnitude of the di¬ 
rection shifts we observe when changing the sky fraction. 

As an independent check, we also produce a cleaned CMB 
map using an internal linear combination (IFC) method. We 
used the HFI maps at 100, 143, and 217 GHz smoothed with 
a 1° FWHM Gaussian kernel. Note that smoothing the data with 
a 1° kernel reduces the solar dipole in the maps by 0.005 %, i.e., 
0.2yuK, which we corrected for afterwards. We then estimate the 
solar dipole amplitude and direction using a Galactic mask that 
removes less than 15 % of the sky to avoid the inner Galactic 
plane where the residuals are most intense. The measurement is 
compatible with the results in Table 2. 

At the end, the amplitude (J) and direction (Galactic longi¬ 
tude, latitude) of the solar dipole measured by Planck-W\ is 

d = 3364.29 + 0.02(stat) + 0.8(sys) + 0.74(FIRAS) yuK 

lon\ _ /263?914 +0?001 (stat) +0?013 (sys) \ 

48?2646 +0?0003(stat) +0?0019(sys )) 

This is to be compared to the official Planck solar dipole 
measurement obtained in combination with the Planck-L¥\'. 

d = 3364.5 + 2.0 yuK 

lon\ _ / 264?00 + 0?03 \ 
lat j " 48?24 + 0?02 ) 

Compared to the WMAP five-year results (J, Ion, lat) = 
(3355 + 8yuK,263?99 + 0?14,48?26 + 0?03; Hinshaw et al. 
2009), this is 9.3 yuK (0.28 %) higher in amplitude while shifted 


by (0.'6,1.'2) in longitude and latitude. Part of the difference 
(0.6 yuK) is due to the revised CMB monopole temperature com¬ 
pared to Mather et al. (1999) (2.725 K). This total dipole (solar, 
orbital, relativistic, and interactions thereof) is removed from the 
calibrated TOI before final mapmaking. 

5.2. Planck-HFI maps 

Frequency maps have been produced using inverse noise weight¬ 
ing. In Figs. 5 and 6 we show the six intensity frequency maps 
from 100 to 857 GHz at full resolution (Aside = 2048). Figure 7 
presents polarization maps at the four first frequencies (100, 143, 
217, and 353 GHz), degraded to lower resolution (Aside = 256) 
in order to enhance the signal-to-noise ratio. Those maps have 
been corrected from bandpass leakage as will be discussed in 
Sect. 7.3.1. In both intensity and polarization, we clearly see the 
emission from the Galactic dust increasing with frequency. In in¬ 
tensity, CMB anisotropies are visible at high latitude in the low- 
frequency channels (between 100 and 217 GHz). In polarization, 
the 100 GHz maps are contaminated in the Galactic plane by 
residual CO leakage coming from bandpass mismatch between 
bolometers. 


5.3. Far sidelobes 

As noted in Planck Collaboration VII (2015), far sidelobes 
(FSFs) affect the response of the instrument to large-scale struc¬ 
ture. In addition, the FSFs also affect the HFI calibration. 

At low frequencies, HFI calibrates by fitting to the si¬ 
nusoidal signal created by the dipole modulated by the 
Planck circular scanning strategy. As outlined in appendix B 
of Planck Collaboration XXXI (2014), this effectively weights 
different parts of the beam in general, and the sidelobes in partic¬ 
ular, by their angle from the spin-axis. For example, far-sidelobe 
contributions close to the spin-axis actually affect the calibra¬ 
tion very little. Similarly, since we are calibrating with signals 
that are “in phase” with the known phase of the main beam as 
it scans, the further a sidelobe contribution is in angle around 
the spin-axis from the main lobe, the less it contributes to the 
calibration. So, a sidelobe contribution that is 90° in scan phase 
from the main lobe, for example, would not contribute to the HFI 
calibration, while something close to the main beam would po¬ 
tentially have a large effect. The change of the gain due to the 
far sidelobes is calculated by fitting the dipole to full timeline 
simulations of the dipole convolved by the FSFs. The factors are 
0.09 % at 100 GHz, 0.05 % at 143 GHz, 0.04 % at 217 GHz, and 
negligible at 353 GHz. The delivered 100-217 GHz maps have 
been scaled by these gain changes. It should be noted that these 
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Fig. 5. Planck-Wl full mission channel intensity maps at 100, 143, and 217 GHz (from top to bottom) after removal of zodiacal 
emission. 
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Fig. 6. Planck-Wl full mission channel intensity maps at 353, 545, and 857 GHz (from top to bottom) after removal of zodiacal 
emission. 
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Fig. 7. Planck-HF\ full mission Q (left) and U (right) polarization maps corrected from bandpass leakage (see Sect. 13.1). from top 
to bottom: 100 GHz, 143 GHz, 217 GHz, and 353 GHz 
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Table 3. Frequency-averaged zodiacal emissivity values for the 
Diffuse Cloud and the three IRAS bands. These are also shown 
in Fig. 8. 


Frequency 

[GHz] 

Cloud 

Band 1 

Band 2 

Band 3 

857 .... 

0.256 + 0.007 

2.06 + 0.19 

0.85 + 0.05 

3.37 + 0.38 

545 .... 

0.167 + 0.002 

1.74 + 0.11 

0.54 + 0.03 

2.54 + 0.18 

353 .... 

0.106 + 0.003 

1.58 + 0.07 

0.39 + 0.02 

1.88 + 0.14 

217 .... 

0.051 + 0.006 

1.30 + 0.07 

0.15 + 0.02 

1.27 + 0.14 

143 .... 

0.022 + 0.010 

1.23 + 0.10 

0.15 + 0.04 

1.16 + 0.22 

100 .... 

0.012 + 0.005 

1.02 + 0.16 

0.08 + 0.03 

0.72 + 0.27 


numbers are uncertain at the 20-30 % level, depending on a mul¬ 
titude of details, such as how the telescope is modelled. 

For the planet photometry, some level of knowledge of the 
amplitude of the FSLs is needed to correctly compare the recon¬ 
structed flux with the planet brightness. However, the relative 
FSL power is lower than 0.3 % (Tauber et al. 2010) for all HFI 
frequencies, which is well below the systematic uncertainties of 
the planet emission models we are using, which are around 5 % 
(see Sect. 3.4). Therefore FSLs can safely be ignored in the 545 
and 857 GHz calibration. 

5 . 4 . Zodiacal emission 

Zodiacal emission is reconstructed and subtracted in 
the same fashion as that used for the 2013 Planck re¬ 
sults (Planck Collaboration XIV 2014). The basic procedure for 
characterizing and removing zodiacal emission from the Planck 
maps is to: 

- make frequency maps for each horn and survey as described 
in previous sections; 

- make survey difference maps for each horn and year; 

- find the date ranges over which each Vside = 256 pixel was 
observed, and veto those pixels that were observed over a 
time-span of more than one week; 

- use the CORE model (Kelsall et al. 1998) to recreate the dif¬ 
ferent zodiacal emission components, assuming blackbody 
emissivities; 

- fit the components to the survey difference maps for each 
horn and year to extract the actual emissivities; 

- use the average of the fitted emissivities to reconstruct the 
implied zodiacal emission seen during each pointing period, 
for each horn, and remove these from each detector. 

The emissivities for each zodiacal component at each of the HFI 
frequencies are given in Table 3 and are plotted in Fig. 8. As 
noted in Planck Collaboration XIV (2014), there seems to be a 
jump between the emissivities for the bands at DIRBE wave¬ 
lengths and the emissivities of Bands 1 and 3 at Planck wave¬ 
lengths. This is being investigated, but is assumed to be a con¬ 
sequence of the assumption in the DIRBE analysis that all three 
bands have the same emissivities, while the Planck analysis al¬ 
lows them to be different. Eor the Planck cosmological studies 
this should be irrelevant, since the zodiacal analysis is being used 
only to remove the interplanetary dust contamination - the over¬ 
all amplitudes of the emissivities, which are completely degener¬ 
ate with the assumed particle density in the bands, are not being 
interpreted physically. 



Fig. 8. Zodiacal emissivities. Data on the left, at wave¬ 
lengths shorter than about 300 yum, are from COBE/DIRBE 
(Kelsall et al. 1998). Data for wavelengths greater than about 
300yum are from Planck (Table 3). In both cases, the blue squares 
represent the emissivity of the Diffuse Cloud. Eor DIRBE, the 
red diamonds represent the fitted emissivity for all three IRAS 
Bands, and the green circles show the values for the Circumsolar 
Ring and Trailing Blob. Eor Planck: the pink, right-pointing tri¬ 
angles are for IRAS Band 3; the brown, left-pointing triangles 
are for IRAS Band 1; and the red, downward-pointing triangles 
are for IRAS Band 2. Eor reference, the lines mark emissivities 
that are unity at wavelengths less than 250yum, but that are pro¬ 
portional to 4“^, /1“\ and 4^ at longer wavelengths. 


6. Noise description and subset differences 

6.1. Map variance 

As demonstrated in Paper 1 , the noise spectra for the Planck-HFl 
bolometers show significant deviation from white noise, result¬ 
ing in correlations between pixels after map projection. At large 
scales, the correlations are dominated by low-frequency noise, 
while at high resolution neighbouring pixels are correlated due 
to time-response deconvolution and filtering. The Planck 2015 
release does not provide a pixel-pixel correlation matrix; only 
the variance per pixel is given for each delivered map. At first 
order, the variance maps are proportional to 1 /Miit, where Ahit is 
the number of samples per pixel. 

The half-difference half-ring maps, projected using the same 
gain, but destriped independently, are a good representation of 
the noise variance in the HEI maps. Indeed, in the difference be¬ 
tween the first and the second half of a ring, the sky signal van¬ 
ishes almost perfectly. Moreover, most of the HEI systematics 
are scan-synchronous and thus also vanish in the difference. 

Table 4 compares the noise per sample from three estimators: 
(a) the mean value of the pre-whitened variance map (i.e., scaled 
using the hit counts); (b) the variance of the pre-whitened half¬ 
ring half-difference map; (c) the average of the half-ring map 
power spectra in the i range 100-5000 (see Sect. 6.2). Eor polar¬ 
ization, the numbers are averages over Q and U for the maps and 
E and B for the spectra. The different estimators are sensitive to 
different kinds of systematic effects, such as time-response resid¬ 
uals and signal gradient in pixels. Nevertheless, the three noise 
estimators give very consistent results. 

The variance maps (/-/, Q-Q, and U-U) are quite inho- 
mogenous, owing to the Planck scanning (which have negligi¬ 
ble wobbling), the relative position of the detectors in the fo- 
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Fig. 9. TT and EE power spectra reconstructed from the half-difference between data subset maps for the dipole-calibrated channels. 


Table 4. Estimation of the noise per sample for intensity (/) 
and polarization {P) estimated from: (a) the variance maps; (b) 
the half-ring difference maps; (c) the pseudo-spectra. Units are 
yuKcMB for 100 to 353 GHz, and MJy sr“^ for the submm chan¬ 
nels. 


Frequency Variance maps Diff. maps Pseudo-spectra 

[GHz] (a) (b) (c) 


1001 . 1538 1531 1410 

loop . 2346 2344 2131 

1431 . 769 758 759 

143P . 1631 1618 1611 

2171 . 1105 1098 1141 

217P . 2512 2486 2440 

3531 . 3692 3459 3780 

353P . 10615 10141 10181 

5451 . 0.612 0.619 0.779 

8571 . 0.660 0.866 0.860 


cal plane, and the rejection of some rings or groups of rings 
(see Paper 1). Moreover, the typical HEALPix pixel size is about 
1.'7 at A^side = 2048 resolution while Planck scans the sky on 
roughly (but not exactly) ecliptic meridians separated by 2.'5 
(near the ecliptic equator). As a consequence, for single sur¬ 
vey maps, lines of empty pixels appear between the scanning 
trajectories around (/,/?) = (0°,+45°). Even when surveys are 
combined, inhomogeneities arise from the HEALPix pixels being 
elongated parallelograms. The axis of their elongation changes 
at the boundaries between the 12 primary HEALPix pixels. In the 
same regions of the sky and in Galactic coordinates, these elon¬ 
gations are parallel with the scanning trajectories, which induces 
moire patterns in the coverage maps. 

The degree of correlation between the Stokes parameters 
within each pixel reflects the distribution of the detector orien¬ 
tations, which results from the scanning strategy. The I-Q and 
I-U correlations are about 14, 9, 6 , and 12% at 100, 143, 217, 
and 353 GHz. The Q-U correlation is about 11, 2, 3, and 8 % at 


100,143, 217, and 353 GHz). In Appendix A.2 we show the sky 
distribution of these correlations. 


6.2. Map differences 

The redundancy of the Planck scanning history and focal plane 
layout provides numerous ways to check data consistency. We 
can create differences between maps built using data splits, as 
described in Sect. 4.3. In the limit that the signal is the same in 
each data subset, the difference map should contain only noise. 
The TOl processing includes several operations that introduce 
correlations on various time scales; these are discussed below. 

In the dipole-calibrated frequency channels (100-353 GHz), 
the signal differences are small enough that the data-split map 
differences can be evaluated at a spectral level, giving insight 
into the residual systematic errors. Eor high-frequency channels, 
we discuss the residuals in the map domain. 

Given maps of two subsets of the data, Ma and Mb, we con¬ 
struct the half-difference as AM = (Ma - Mb) /2. We compute 
the power spectrum in temperature and polarization of the half¬ 
difference, masking the sky with the Planck point source mask 
as well as the galaxy masks used in Planck Collaboration XI 
(2015), i.e., leaving 65, 59, 48, and 32% of the sky unmasked 
at 100, 143, 217, and 353 GHz, respectively. 

In order to use this half-difference map to assess noise in 
the full maps, we account for widely varying integration time 
in the two subsets using a pixel-by-pixel weight map, which 
is multiplied by the half-difference map AM prior to comput¬ 
ing the a ngular power spectrum. The weight is constructed as 
W = 21 where ^a is the hit count map 

for Ma and ^5 is the hit count map for M 5 . In the limit that the 
half-difference map consists entirely of white noise, this exactly 
accounts for the differences in the hit counts. The TT and EE 
spectra of the difference maps are plotted in Eig. 9 and are de¬ 
scribed in the sections below. The BB spectra are nearly the same 
as the EE spectra and are not shown. 
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6.2.1. Half-ring map differences 

The half-ring difference is sensitive to high-frequency noise, 
since most low-frequency modes (on time scales longer than 1 h) 
are common to both data sets and thus vanish. In the harmonic 
domain, the noise is nearly white with an amplitude compati¬ 
ble with the noise estimated in the map domain (see Sect. 6.1). 
At large multipoles, the noise blows up due to the time trans¬ 
fer function deconvolution, before being cut off by the low- 
pass filter. At lower multipoles, half-ring differences show low- 
frequency noise residuals due to the destriping. Indeed, the de- 
striping is performed independently for each half, or essentially 
half the data are used to solve the offsets for the full ring maps. 
The residuals from the offset determination are therefore ex¬ 
pected to be twice as large as in the full-mission map. 

In addition, the deglitching operation performed during the 
TOI processing uses the full data set to estimate the signal in 
each ring, thereby introducing some correlation between the two 
halves of each ring. Taking the difference between the two half¬ 
rings in fact removes the correlated portion of the noise at the 
few percent level. 

6.2.2. Half-mission map differences 

With half-mission differences, we can check for long-time-scale 
variations and for apparent gain variation with time due to ADC 
nonlinearities. Moreover, due to slightly shifted pointing be¬ 
tween the first and second halves of the mission, the effect of 
a signal gradient within a pixel (especially on the Galactic plane 
where the signal is strong) is larger than for the half-ring map 
differences. 

Because the number of observations in a given pixel can 
be very different between the two half-mission maps, using a 
weighting as described above is essential. Including the weight¬ 
ing, the half-mission differences give a power that is 10-20% 
higher than the corresponding half-ring difference. This fraction 
of additional power is nearly the same in all the channels 100- 
353 GHz, and is the same in both temperature and polarization. 
We understand this small additional power to be due to effects 
from the TOI processing that introduce correlations in the noise 
between the subsets. The half-ring maps, as stated above, have 
correlations introduced by the deglitcher that are subsequently 
removed by the differencing. These correlations are not present 
between the half-mission data sets, so their difference shows a 
higher noise power. 

6.2.3. Detector-set map differences 

This difference probes systematic effects that are bolometer- 
dependent. Note that, in the case of 143, 217, and 353 GHz, 
the detector-set split excludes the unpolarized detectors, and the 
noise in TT is correspondingly higher than in the half-mission 
and half-ring split. The 100 GHz channel has only polarization- 
sensitive bolometers and the TT spectrum of the difference is 
much closer to the spectrum seen for the other data splits. 

There are several other effects that make the power spectrum 
of the detector-set difference stand apart. A unique time response 
function is deconvolved from each bolometer. In the half-ring 
and half-mission data splits, the deconvolved function is identi¬ 
cal between the two halves. However with the detector-sets, the 
time response is in general slightly different between the two 
halves. This effect leads to a tilt in the spectrum of the detector- 
set difference maps relative to the half-ring or half-mission split. 



10 ° 10 ^ 10 ^ 


Fig. 10. EE power spectra reconstructed from the half-difference 
between data subset maps for the dipole-calibrated channels at 
low multipoles compared to the noise estimation from the FFP8 
simulations. CMB signal from Planck 2015 is plotted in dashed 
lines. 


Moreover, at 353 GHz, signal residuals are larger due to relative 
calibration uncertainties between detectors. 

6.2.4. Map differences at low-^ 

At low multipoles, despite the huge progress in the control of the 
systematics, data are still contaminated by systematic residuals. 
Figure 10 shows the EE power spectra from the half-difference 
maps at 100, 143, and 217 GHz and compared to the noise power 
spectrum from FPP8 simulations. The half-ring differences are 
compatible with noise while, at multipoles typically lower than 
50, detector-set and half-mission differences are dominated by 
excess power which is larger than the EE CMB signal. The ori¬ 
gin of the excess power will be explored in a forthcoming publi¬ 
cation. 

6.2.5. High frequency channels 

For the highest frequency channels (545 and 857 GHz) the data- 
split map differences are dominated by residual signal. Figure 11 
shows the rms of the differences of intensity maps at 545 and 
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Fig. 11. Rms of the residual signal in difference maps at 545 and 
857 GHz, as a function of signal level in the full map. The solid 
coloured curves show the rms of the data, while dashed coloured 
curves show the rms of a simulated noise map. The diagonal 
dotted lines indicate 1 % and 10 % of the signal. 


857 GHz for half-ring, half-mission, and year data splits com¬ 
pared to the same data split performed on a simulated noise 
map. At low signal, the difference is consistent with instrumental 
noise. At high signal levels, an additional residual appears in the 
difference map that is roughly proportional to the signal level. 
Part of this is due to pointing errors. For year and half-mission, 
the effect is enhanced by the combination of residual gain vari¬ 
ations and the relative difference of pointing between the two 
splits. Over most of the sky, the signal is reproducible to better 
than 1 % for these frequencies. Bolometer map differences (not 
shown here) are, in addition, sensitive to the relative calibration 
error. 

6.3. Noise cross-correlation 

Here we check for correlations in the noise by computing cross¬ 
spectra between the difference maps described earlier. We look 
at 100 (Fig. 12), 143 (Fig. 13), and 217 GHz (Fig. 14) in com¬ 
parison with the expectations from projecting noise realiza¬ 
tions on the sky (using the FFP8 noise realizations described 
in Planck Collaboration XII 2015 and the end-to-end simula¬ 
tions described in Paper 1). 

When the half-mission cross-spectra of half-ring differences 
are computed, the results are roughly consistent with the FFP8 


noise simulations. At 143 GHz in temperature, the end-to-end 
simulation produces a slight rise in power at low multipoles that 
is not seen in the data. 

Large correlations are seen in the half-ring cross-spectra of 
half-mission differences. These are at least partially induced by 
our processing since the end-to-end simulations also show cor¬ 
relations that are not as large in amplitude as those seen in the 
data, but show a similar spectral shape. These correlations are 
mainly due to the deglitcher, as described above. 

6.4. Temperature-polarization cross-variance 

In absence of spatial correlations, noise correlations between 
temperature I and polarized Q and U modes vanish in the 
harmonic domain, thanks to the orthogonality of the spherical 
harmonic decomposition. Consequently the TE and TB auto¬ 
spectra are not biased by noise in the way that the TT, EE, and 
BB spectra are. In practice, transfer function deconvolution, fil¬ 
tering, and pixelization effects can produce spatial correlations at 
high multipoles, resulting in a noise bias that is observed in the 
TE and TB angular power spectra. In Fig. 15 we compare the 
auto and cross-spectra for the half-ring, half-mission, and detec¬ 
tor set splits. These pseudo-spectra have been built by masking 
Galactic emission and point sources (approximately 40 % of the 
sky). The auto-spectra are biased at high multipoles (starting at 
i ^ 1500). The amplitude of this bias and its sign depend on 
the frequency and on the mode considered. Nevertheless, none 
of the cross-spectra show significant departures from the null ex¬ 
pectation. 

We observe that the amplitude of the noise bias in auto¬ 
spectra is mitigated when adding more independent data sets, 
such as detectors or surveys (survey maps show larger ampli¬ 
tude than half-mission and full-mission). These results are fully 
reproduced in the FFP8 simulations. 

6.5. Subset map cross-covariance 

When mapping subsets of the available data (selecting detec¬ 
tors and/or time spans) we have a choice between solving for 
independent baseline offsets for the subset in question or reusing 
full-mission, full-frequency baselines (as in the 2015 HFI map 
release). Full-mission baselines are more accurate, leaving less 
large-scale noise in the maps, but introduce noise correlation be¬ 
tween detector-set maps. 

We can measure the resulting bias in cross-spectra through 
noise simulations. Comparison between noise spectra from a 
Monte Carlo analysis at 100 GHz using both methods of destrip- 
ing is shown in Fig. 16 for the case of detector sets. Noise spectra 
for each detector set show more large-scale power when using 
independent baselines than in the full mission case. On the other 
hand, the cross-power spectrum is biased by up to a few times 
10“^yuK^ below ^ = 10; the same is true for the half-mission 
subset. 


7. Systematic effects 

We now describe the major systematic effects that could poten¬ 
tially affect the maps: the gain variations; errors in the abso¬ 
lute gain determination; errors in the polarization efficiency and 
orientation; and, most of all, the detector-to-detector gain mis¬ 
match. The latter includes bandpass mismatch (which affects the 
response to foregrounds of the detectors at the same frequency) 
and relative gain uncertainties, both of which create intensity-to- 
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Fig. 12. 100 GHz difference map cross spectra. Left: half-mission (HM) correlation of half-ring differences (HR). Right: half-ring 
(HR) correlation of half-mission difference (HM). The real data are red dots. The end-to-end simulation are black stars. One and 
two sigma contours from ten FFP8 noise realizations are shaded grey. 
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Fig. 13. Same as Fig. 12 for 143 GHz difference map cross spectra. 


polarization leakage. All these effects are constrained using tests 
involving the combination of maps, residuals in maps, cross¬ 
power spectra, and dedicated simulations. 


7.1. Gain stability 

Gain stability has been significantly improved with respect to 
the Planck 2013 release. This is mainly due to the ADC correc¬ 
tion, combined with the new determination of the time transfer 
function (see Paper 1). The amplitude of the apparent gain vari¬ 
ation has been improved from 1-2 % to less than 0.5 % for all 
cases. Residual gain variations are compatible with zero when 
including the correction for the long-time-constant residuals, as 
discussed in Sect. 3.3. 


We check the stability of the gain over time using the same 
tool as in Planck Collaboration VIII (2014), called bogopix. For 
each bolometer, the code fits simultaneously for the gain gr and 
the offsets Oy for each ring, marginalizing over the sky signal T: 

Pt — §r ^ iA^tp ■ Pp F ^orb) + • (15) 

Given the low amplitude of the observed gain variations (less 
than 0.5 %), we linearize Eq. (15) and solve by iteration (see 
Planck Collaboration VIII 2014); one or two iterations are suffi¬ 
cient to ensure convergence. To initialize the iterations, we start 
from the constant gain solution G (see Sect. 3.2). 

We compute the gain variations from single-bolometer data 
(neglecting polarization). Polarized signals will affect the gain 
determination. To reduce this bias, we ignore sky regions where 
the polarized emission is the strongest, which lie mostly in the 
Galactic plane. 
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Fig. 14. Same as Fig. 12 for 217 GHz difference map cross spectra. 


Figure 17 shows the results of bogopix for bolometers at 
100, 143, and 217 GHz, smoothed over a 4-day period. At higher 
frequencies (353 GHz and above), the gain variations are much 
lower than the gain uncertainty. Owing to the Planck scanning 
strategy, the Galactic foreground is larger for some rings, whilst 
the orbital dipole amplitude is almost constant. This increases 
the dispersion around those regions and potentially induces some 
bias in the gain determination. In the end, we find gain variations 
with amplitudes lower than 0.3, 0.4, and 0.5 % at 100, 143, and 
217 GHz, respectively. 

The residual apparent gain variations are essentially com¬ 
ing from: the uncertainty in the current ADC correction; the un¬ 
certainty in the long-time-constant estimation; and the effect of 
long-term thermal variations on the bolometer and electronics 
response. 

7.2. Calibration accuracy 
7.2.1. Inter-frequency accuracy 

The precision of the calibration can be assessed by looking for 
residual dipoles in the maps. If the calibration of a map is slightly 
incorrect, the removal of the solar dipole in the mapmaking pro¬ 
cess leaves a residual dipole. However, identifying such a resid¬ 
ual dipole is difficult because of the presence of other sources of 
power at 7” = 1, mostly due to Galactic emission and zodiacal 
light, but also to imperfect correction of systematic effects such 
as far sidelobes. 

Following the method presented in Sect. 5.1 used to esti¬ 
mate the solar dipole direction and amplitude, we cleaned the 
Galactic emission from the HFI maps using a local correlation 
with the 857 GHz map. Adding the solar dipole that was re¬ 
moved in the mapmaking process (Sect. 5.1) to Dj-es produces 
a map that contains the true solar dipole. We then fit for its 
amplitude, fixing its direction to the official Planck value (Ion, 
lat = 264?00,48?24) to limit the effect of other residuals still 
present in Dj-es (Galactic, systematic effects). The fit is done us¬ 
ing sky pixels where hsi < 2 (faintest 37 % of the sky) to limit 
the effect of Galactic emission residuals (CO, free-free emission, 
and small scale dust SED variations). Table 5 gives the ratio of 
fitted amplitude to the removed dipole found for the HFI maps 


Table 5. Ratio of amplitudes of the fitted dipole (A^t) and the 
removed dipole (Ann = 3364.5 jjK). The direction of the dipole 
removed from the data (lon=264?00, lat=48?24) was constrained 
here. The fit was computed for different sky fractions from 30 to 
70 %. Statistical and systematic uncertainties are also indicated. 


Frequency Afit/Aj-m Statitiscal Systematic 
[GHz] uncertainties uncertainties 


100 . 1.00010 +0.00006 +0.0001 

143 . 0.99988 +0.00012 +0.0001 

217 . 1.00184 +0.00027 +0.0003 

353 . 1.00568 +0.00185 +0.0020 

545 . 1.02515 +0.01627 +0.0190 


at frequencies from 100 to 545 GHz. Because effects other than 
a miscalibration can contribute to a residual dipole in the maps, 
these ratios provide upper limits on the calibration accuracy at 
each frequency. The results indicate that the calibration at 100 
and 143 GHz is precise at a level of few lO""^. At 217 GHz, the fit 
is compatible with a residual dipole at the 0.2 % level. At higher 
frequencies, the fits indicates residuals at 0.52 % and 1.23 % at 
353 and 545 GHz, respectively. 

These results are in agreement with those obtained while 
performing component separation, as shown in Table 4 from 
Planck Collaboration IX (2015) and in Planck Collaboration X 
(2015). They are also in agreement with the results from the cos¬ 
mological parameter determination, where intercalibration coef¬ 
ficients are also fitted for (see Planck Collaboration XI 2015). 
The agreement between these measurements computed over dif¬ 
ferent multipole ranges highlights the quality of the Planck-Wl 
calibration, together with the accuracy of the transfer function 
reconstruction. 

7.2.2. Intra-frequency accuracy 

For polarization reconstruction with Planck-Wl data, we have 
to combine data from several detectors. Any relative calibra¬ 
tion error will induce an intensity-to-polarization leakage (see 
Sect. 7.3). For the CMB channels, we have assessed the relative 
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Fig. 15. Pseudo-power spectra for TE (/eft) and TB (right) 
for each frequency (from top to bottom: 100, 143, 217, and 
353 GHz). The auto-spectra are shown in black. Cross-spectra of 
half-ring (HR), half-mission (HM), and detector-set (DS) half¬ 
differences are shown in blue, red, and green, respectively. A 
Galaxy and point source mask, leaving 40 % of the sky, was used 
in all cases. 



Fig. 16. Average EE spectra obtained from 100 simulations of 
detector-set noise maps. The latter have been produced using ei¬ 
ther the full-mission baseline destriping (MBD) or the destrip- 
ing run on each subset independently (run baseline destriping or 
RBD). Both the increased low-^ noise in the RBD case and the 
increased cross-spectrum noise bias in the MBD case are appar¬ 
ent. 


calibration accuracy for each detector at a given frequency using 
two complementary methods. 

As in section 6.2 (see figure 14) of Planck Collaboration VIII 
(2014), we derive relative inter-calibration factors for each de¬ 
tector (for 100 to 353 GHz), rescaling their cross-pseudo-power 
spectra, estimated over 40% of the sky (30% for 353 GHz) in 
the i range 25-300, which encompasses the first acoustic peak. 
We used colour-correction factors at 353 GHz, because even at 
high latitude the dust emission is large. As in 2013, we keep the 
maximum of these factors as a conservative estimate of the rela¬ 
tive calibration accuracy. In 2015, we find 0.09, 0.07, 0.16, and 
0.78 % for 100, 143, 217, and 353 GHz, respectively (compared 
to 0.39, 0.28, 0.21, and 1.35% in 2013). Since single-detector 
maps are built ignoring polarization, these values should be con¬ 
sidered as conservative upper limits on the relative detector-to- 
detector calibration accuracy. 

We complemented these estimations by analysing the solar 
dipole residual on the differences of single detector maps. We 
fit the dipole amplitude fixing its direction while masking 30 % 
of the sky in the Galactic plane to avoid regions affected by 
band-pass differences. We find maximal amplitudes 0.5, 0.6, and 
3.0 yuK for 100, 143, and 217 GHz, respectively, which - relative 
to the solar dipole amplitude (3364.5yuK, see Sect. 5.1) - gives 
accuracies of the same order as the aforementioned spectra anal¬ 
ysis. As in this previous method, the main limitation comes from 
polarization which is ignored in the single detector maps. 

While significantly better than for the 2013 release, calibra¬ 
tion mismatch between bolometers at a given frequency is one 
of the main systematic residuals contaminating the HFI large an¬ 
gular scales in polarization, as explained in the next section. 

7.3. Intensity-to-polarization leakage 

Any gain mismatch between the measurements of detectors be¬ 
longing to the same frequency channel will result in intensity-to- 
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Fig. 17. Gain variation with ring number for each bolometer 
estimated using bogopix. From top to bottom: 100, 143, and 
217 GHz. Gain values for individual rings (gray dots) have been 
smoothed with a 4-day width. The gain variations are lower than 
0.3, 0.4, and 0.5 % at 100, 143, and 217 GHz respectively. 


polarization and cross-polarization leakage in the channel maps. 
In Planck, the dominant leakage effect has three main origins: 

- monopole mismatch from the uncertainty in the mean offset 
determination; 

“ gain mismatch that produces leakage from the whole inten¬ 
sity signal into polarization; 

- bandpass mismatch (hereafter BPM) that mainly generates 
intensity-to-polarization leakage from foreground emission 
(with a non-CMB spectrum). In the case of HFI, the leakage 
effect is dominated by CO and thermal dust emission. 

All these leakage sources are especially important for the large 
angular scales. Beam-mismatch polarization leakage occurs at 
small angular scales and is discussed in Paper 1. Although the 
first two mismatches can be minimized by obtaining more accu¬ 
rate measurements of offsets and gains respectively, the BPM 
cannot be removed in the mapmaking process if we want to 
project CMB and foregrounds at the same time. 


The power absorbed by a given bolometer b at time t is ex¬ 
pressed using the Stokes parameters {Ip ,Qp, Up) which char¬ 
acterize the emission in intensity and polarization in the corre¬ 
sponding sky pixel p. The polarized HFI channels are calibrated 
using the CMB orbital dipole and the total calibrated power ab¬ 
sorbed by the bolometer b can be written as 

= (1+i.n) 

X Ci(l + 4p,t) + P* {Q\ cos (fi + ul sin <^f)] 

+ Coffset + «*’ (16) 

where the polarization efficiency and the polarization angle 
= 2{il/t + a^) are explicitly dependent on the bolometer b, 
with the index k ranging over the different sky components. 
Additionnaly we have the following definitions: 

- encodes the gain mismatch of bolometer b with respect 
to the mean calibration of the channel; 

- corresponds to the overall offset of bolometer b, which 
is small but not vanishing; 

- Ck is the average transmission of sky component ^ in a 
given channel and is the bandpass mismatch specific 
to bolometer b, and affecting all sky components except the 
CMB (see description below). 

Each of these e terms is responsible for leakage from intensity 
to polarization in a manner that can in principle be quantified 
and corrected for, as described hereafter. Note that we only con¬ 
sider first-order terms in 6, as any higher-order contribution is 
negligible. 


7.3.1. Bandpass mismatch (BPM) 

Each emission component k (where k = CMB, dust, synchrotron, 
etc.) is integrated over the bandpass of the detector according to 
a given spectrum /^(v). Since the polarized HEI channels are cal¬ 
ibrated using the CMB orbital dipole, we define the transmission 
coefficients 

f fcMB(y)^v^y 
= Cjc(1+4p,/c)’ 

where is the spectral response of bolometer b and Q = 
Y,b C'^/^boio is the average value of the in a given channel. 
These transmission coefficients express the ^-component emis¬ 
sion in CMB units. If all bolometers had the same spectral re¬ 
sponses then 6gp^ would be equal to zero, i.e. = Q, in which 
case no BPM-related leakage would be produced. 

Considering only bandpass mismatch corresponds to setting 
^gain “ ^ ^offset “ ^ in Eq. (16). Then, ordering all the data 
samples, mf, for a bolometer observing a position p on the sky 
into a single vector defining A to be the pointing matrix in 
temperature and polarization, and n the noise vector, Eq. (16) 
reads 



Using all bolometers b within a channel, the mapmaking proce¬ 
dure solves for the total signal Stokes parameters (/^ ^, 2^ ^0 
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in pixel p, formally computing, 
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( 20 ) 


wherer^ = (A^N A^N A. We have introduced the ma¬ 
trix A^, the elements of which are equal to zero except for the 
diagonal elements relevant to bolometer b, which are set to 1. 
The last term of Eq. (20) is the leakage term in pixel p, where 
intensity will leak into Q and U, Q into 1 and U, and U into / and 
2, according to the mismatch coefficients 6gp ^ and the values of 
the 3x3 matrix 
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Considering all pixels, the quantities correspond to nine 
sky maps for bolometer b. These maps can be fully determined 
from the mapmaking solution and may be understood as patterns 
of the mismatch leakage. In practice, cross-polarization leakage 
and polarization-to-intensity leakage are negligible compared to 
the intensity-to-polarization contribution and we therefore con¬ 
sider the latter only. The T jq and T lu maps have been system¬ 
atically produced by the mapmaking pipeline.^ With these as¬ 
sumptions, the BPM-induced leakage in Q and U for the sky 
component k reads 


rBP,^ 

^IQJU 


= CkP 


^bolo 1 

/ V ^BP,^ ^ IQ,IU • 
b=0 


( 22 ) 


In consequence, for a given calibrated intensity template of 
the sky component k (i.e., = Qf^) we may compute leak¬ 

age correction maps as 


^bolo 1 ^bolo 1 

jCon:,k _ jk \ ^ b pZ? _ jk \ ^ ^ pZ? 

^IQ,IU ~ ^template /_j ^BP,^ ^ IQ,IU “ ^template Z_^r ’ 

b=t) b=0 

(23) 

where the last equality uses the fact that Yjh=o ^ /t/ “ ^ 

construction. 

Leakage correction maps have been produced for all po¬ 
larized HFI channels. The relevant foregrounds at these fre¬ 
quencies are dust (all channels) and CO (all channels except 
143 GHz). To do so, the coefficients have been computed 
from Eq. (17), where the spectral responses of the bolometers 
H^(y) are those obtained from pre-launch ground-based mea¬ 
surements of the bandpasses (Planck Collaboration IX 2014). 
The dust spectrum is taken as a greybody with spectral index 
P = 1.62 and temperature T = 19.7 K, which are the all-sky av¬ 
erage values found in Planck Collaboration XI (2014). For the 
intensity template required in Eq. (23), we use the thermal dust 

^ The TII pattern map quantifies the correction that should, in prin¬ 
ciple, be brought to the I channel map, given that some intensity has 
leaked into polarization. The correction is, however, negligible and not 
taken into account here. The same is true for Tqq and Tuu- 


intensity maps at 353 GHz obtained from the Planck thermal 
dust model (Planck Collaboration XI2014). Combining all these 
ingredients, the dust correction maps are produced ac¬ 

cording to Eq. (23) and delivered in the 2015 HFI data release 
(Fig. 18). 





Fig. 18. Dust leakage correction maps from ground-based mea¬ 
surements of the bandpass in Q (left) and U (right) at all HFI 
channels: 100, 143, 217, and 353 GHz (from top to bottom). 


Note, however, that the reliability of these corrections is lim¬ 
ited by uncertainties both in the physical nature of the fore¬ 
ground components and in the determination of the bolometer 
spectral responses. For the sake of simplicity, a constant spec¬ 
tral index and a constant temperature across the sky have been 
assumed for the thermal dust emission. Furthermore, the cali¬ 
brated thermal dust intensity templates are those derived from 
the 2013 Planck thermal dust model which, while close, does not 
strictly correspond to the calibration of the 2015 maps. Also, the 
leakage corrections are particularly sensitive to the differences 
in transmission between bolometers (i.e., the coefficients); 
small uncertainties on those will yield large uncertainties in the 
final correction maps. In conclusion, the bandpass leakage cor¬ 
rections should not be taken at face value, but should be thought 
of as order-of-magnitude estimates only. We only advocate the 
use of these correction maps to test the stability and estimate un¬ 
certainties of any further results using the HFI polarization maps. 
A result solely obtained by applying the corrections will not be 
reliable. 

7.3.2. Calibration and monopole mismatches 

Using the same formalism as above, calibration mismatch is 
computed by setting = 0 and = 0 in Eq. (16). 
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Following closely Sect. 7.3.1, one finds that the total intensity- 
to-polarization leakage due to calibration mismatch is 

^bolo~l 

^1q,IU - Xi ^ ^ 4ain 

k b=0 

^bolo~l 

- /dipole X (24) 

b=0 

where at first order, for low-frequency maps, the solar dipole sig¬ 
nal (k = dipole, Cdipoie = 1 by construction) provides the domi¬ 
nant contribution to the calibration mismatch leakage effect. 

Setting 6gp^ = 0 and = 0 in Eq. (16), one shows in a 
similar fashion that the monopole intensity-to-polarization leak¬ 
age is simply 


b^holo 1 

rmono _ b yb /orx 

^IQ,IU - ^offset ^ IQ,IU ’ 

b=0 

where the monopole mismatch is modelled using a constant sky 
template while the amplitude of the mismatch is 

encoded in 

Although the BPM coefficients 6gp can be evaluated di¬ 
rectly from foreground modelling (assuming a given spectrum 
of the dust and using the spectral responses of the detectors), 
this is not the case for 6^^.^ and It is therefore not possible 
to provide correction maps for these leakage effects by comput¬ 
ing Eq. (24) and (25) directly. However, one may consider the 
possibility of fitting these quantities from the maps themselves, 
by using the x E^^ and as templates of the gain 

and monopole leakages respectively. Such a method, dubbed the 
“generalized global fit” (GGE), has been implemented and is fur¬ 
ther described in Appendix A.3.2. 

The leakage maps produced with the GGE method for the 
353 GHz channel are delivered in the 2015 release and corrected 
for BPM, calibration, and monopole leakage simultaneously. 


- we convolved these maps with the single bolometer, single 
survey effective beams using EEBeCoP (Mitra et al. 201 1); 

- using the Crab IRAM map as a template, and the polarization 
angles a, we modelled the intensity map described above in 
the Crab region, as a function of an angular offset Aa. 

We then fitted for the values of the angular offset Aa. To do 
that, we first removed the background from the single bolometer 
maps. We built a noise model combining the single detector pixel 
variance with the noise of the IRAM observation, taking into 
account the smoothing applied. We used the Rosset et al. (2010) 
values as a prior. The resulting angular offsets are presented in 
Pig. 19. Corrections are compatible with zero, and this analysis 
does not favour an update of the ground-based parameters. 
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Fig. 19. Estimated angular offset from comparison of the sin¬ 
gle bolometer, single detector maps with the IRAM Crab nebula 
maps, combining the first four surveys. 


7.4. In-flight validation of the polarinneter efficiency and 
orientation 

As discussed in Sect. 4.5, the polarimeter efficiency and ori¬ 
entation used in this release are taken from ground measure¬ 
ments (Rosset et al. 2010). In order to validate these num¬ 
bers in flight, we used the Crab nebula maps obtained with 
the IRAM 30 m telescope and the 90 GHz XPOL polarimeter 
(Aumont et al. 2010). These maps consist of /, Q, and U mea¬ 
surements with an angular resolution of 27" of a 10'-wide region 
around the Crab nebula (Tau A, Ml, or NGC 1952, at J2000 
coordinates RA = 5^34”^32^ and Dec = 22°00'52"). The same 
region was observed by Planck once per survey, with different 
scan directions for odd and even surveys. We compared single 
survey, single bolometer maps of the Crab region with a model 
obtained from the IRAM maps, and solved for the best values of 
polarimeter angle and efficiency. 

Prom single survey, single bolometer data, we can only solve 
for an intensity map, which projects on the sky the total power 
Pt described in Eq. (1). This power depends on the true value of 
the polarization angle a specific to the detector. 

We compared the single bolometer, single survey maps with 
a model obtained with the following procedure: 

- we pixelized the IRAM observations on a HEALPix grid with 
A^side = 2048, rotating to Galactic coordinates; 


We used the same procedure also to fit the polarization effi¬ 
ciency p, but the result is completely dominated by the ground- 
based calibration prior. 

If we assume that the CMB anisotropies have vanishing TB 
and EB power spectra, i.e., that there are no parity violating 
physical mechanisms in the early Universe, we can also check 
whether the overall polarizer angle of Planck-Y{F1 is compati¬ 
ble with zero. Planck Collaboration XXXI (2015) show that the 
CMB TB and EB spectra measured by HEI are consistent with 
zero. Their analysis gives a polarizer angle within 0.3° of zero, 
which is identical to the systematic error of the ground-based 
measurements. This is a factor of five improvement over the 
WMAP final results (Hinshaw et al. 2013) and comparable with 
ACT (Naess et al. 2014; see also the review by Gubitosi & Paci 
2013). 

8. Conclusions 

This paper has described the processing applied to construct the 
Planck-HFl maps delivered in the 2015 release. It has also as¬ 
sessed the main characteristics of the maps in terms of noise and 
systematics, in particular resulting from ADC corrections and 
bolometer long time constants. Since the last release, the calibra¬ 
tion has been upgraded and is now significantly more accurate. 
At low frequency, it is now independent and based on the orbital 
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Table 6. Main characteristics of HFI Full Mission Maps. 


Quantity 







Notes 

Reference frequency v [ GHz]. 

100 

143 

217 

353 

545 

857 

al 

Number of bolometers. 

8 

11 

12 

12 

3 

4 

al 

Effective beam solid angle Q [arcmin^]. 

106.22 

60.44 

28.57 

27.69 

26.44 

24.37 

bl 

Error in solid angle ctq [arcmin^] . 

0.14 

0.04 

0.04 

0.02 

0.02 

0.02 

bl 

Spatial variation (rms) AQ [arcmin^]. 

0.20 

0.20 

0.19 

0.20 

0.21 

0.12 

b3 

Effective beam EWHMi [arcmin] . 

9.68 

7.30 

5.02 

4.94 

4.83 

4.64 

b4 

Effective beam FWHM 2 [arcmin] . 

9.66 

7.22 

4.90 

4.92 

4.67 

4.22 

b5 

Effective beam ellipticity e . 

1.186 

1.040 

1.169 

1.166 

1.137 

1.336 

be 

Variation (rms) of the ellipticity Ae . 

0.024 

0.009 

0.029 

0.039 

0.061 

0.125 

bl 


Sensitivity per beam solid angle fuKcMs] . 7.5 4.3 8.7 29.7 

[kJy sr"^]. 

Temperature Sensitivity fuKcMB deg]. 1.29 0.55 0.78 2.56 

[kJy sr"^ deg]. 

Polarization Sensitivity fuKcMB deg]. 1.96 1.17 1.75 7.31 


Calibration accuracy [%] . 

0.09 

0.07 

0.16 

0.78 

l.l(+5) 

1.4(+5) 

d 

CIB monopole prediction [ MJy sr“^]. 

0.0030 

0.0079 

0.033 

0.13 

0.35 

0.64 

el 

Zodiacal light level correction [Kcmb]. 

[MJysr'] _ 

4.3 X 10-’ 

9.4 X 10-’ 

3.8 X 10-« 

3.4 X 10-5 

0.04 

0.12 

el 

el 


Channel map reference frequency, and channel identifier. 

Number of bolometers whose data were used in producing the channel map. 

Mean value over bolometers at the same frequency. See Sect. 4.2 in paper 1. 

As given by simulations. 

Variation (rms) of the solid angle across the sky. 

FWHM of the Gaussian whose solid angle is equivalent to that of the effective beams. 

Mean FWHM of the elliptical Gaussian fit. 

Ratio of the major to minor axis of the best-fit Gaussian averaged over the full sky. 

Variability (rms) on the sky. 

Estimate of the noise per beam solid angle as given in bl . 

Estimate of the noise in intensity scaled to 1° assuming that the noise is white. 

Estimate of the noise in polarization scaled to 1° assuming that the noise is white. 

^ Calibration accuracy (at 545 and 857 GHz: the 5% accounts for the model uncertainty). 

According to the Bethermin et al. (2012) model, whose uncertainty is estimated to be at the 20 % level (also for constant v/y). 

Zero-level correction to be applied on Zodical light corrected maps. 

dipole signal, while the planets Uranus and Neptune are used to 
calibrate the high end of HFI, achieving 6.1 and 6.4 % absolute 
photometric calibration at 545 and 857 GHz, respectively. This 
has allowed us to measure a consistent CMB solar dipole with an 
unprecedented accuracy better than 10“^ and in agreement with 
the independent determination by LFI. 

Table 6 gives a quantitative assessment of the main charac¬ 
teristics of the Planck HFI maps from the 2015 release. They 
now cover the entire Planck HFI cold mission (885 days). The 
HFI aggregated sensitivity (referring to a weighted average of 
the 100, 143, and 217 GHz channel maps) is 26yuKcMB-arcmin 
in temperature and 52 yuKcMB-arcmin in polarization. 

The noise in the maps shows some small low-frequency ex¬ 
cess on top of white noise prior to time constant deconvolution. 

The latter then naturally raises the higher part of the noise spec¬ 
tra in the multipole domain. We have identified a low-level noise 
correlation in particular between half-ring and detector subsets 
that is not directly reproduced by simulations, although the level 
is small compared to the CMB signal. 

The raw sensitivity must be matched by a long list of con¬ 
straints on any possible systematic effects. This list includes: an 
absolute calibration at a level of 0.1 % to 1.4 % depending on the 
frequency; a resulting apparent gain variation of less than 0.5 %; 


and a knowledge of the polarization angle and polarization ab¬ 
solute value respectively, at the degree level and the 1 % level. 
The instrumental beam has been measured at the percent level 
by using multiple planet crossings. 

Despite the huge progress made in the understanding of 
all the aforementioned systematic effects, Planck-R¥\ polariza¬ 
tion maps are still dominated by systematic residuals at large 
scales. These are essentially coming from the temperature-to- 
polarization leakage resulting from the mismatch between the 
bolometers that are combined to reconstruct linear polarization 
maps. The origins of the leakage effects include: mismatch of the 
zero level from uncertainty in the offset determination; mismatch 
from gain uncertainty (even at the 10“^ level); and bandpass 
mismatch. Corresponding first order corrections for monopole, 
dipole, and bandpass mismatch are provided (as described in 
Appendix A.3. 2) but residuals are still found to be larger than 
noise at very large scales. As a consequence, the Planck-RPl 
polarization maps at large scales cannot yet be directly used for 
cosmological studies. 
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Appendix A: HFI map product description 

Here we summarize the HFI map products that are part of the 
Planck 2015 data release. 

A.1. Map products 

The 2015 release contains many different maps whose details 
are described in the subsections below. All maps are given as 
HEALPix vectors with NESTED ordering, in Galactic coordi¬ 
nates, at the resolution corresponding to Aside = 2048 (for high 
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resolution maps). Depending on the type of product, these vec¬ 
tors are packaged into a binary table and written into a FITS file. 
The table contains in most cases 50 331 648 rows (the length of 
the HEALPix vector for A^side = 2048) and either three columns 
(/, //, H, respectively for intensity, noise (co)variance, and hits) 
for the temperature-only cases or 10 columns (/, Q and U signals 
plus the six various /, 2, and U noise (co)variances and hits) for 
the polarized cases. Pixels with a condition number larger than 
10^ see their hit number brought down to zero. A pixel with zero 
hit has an intensity value of -1.6375 x 10^^ (the HEALPix con¬ 
ventional null value). 

The general matrix of products is pictured in Fig. 4. The 
main products are the maps of the full channels, covering the full 
mission. For characterization and analysis purposes, the chan¬ 
nels have been split into independent detector sets. The detec¬ 
tor sets for the polarized channels are groups of four PSBs that 
can be used to build a sky map in temperature and polarization. 
Bolometers insensitive to polarization (SWB) of the same fre¬ 
quency channels are not used in the detector sets. The detector 
set names and the bolometers they use are listed in Table A. 1 . 
For completeness, full mission maps built using each of the un¬ 
polarized bolometers are also provided. 

The mission duration has been split into single surveys, 
“years” and the two halves of the mission duration. The sur¬ 
veys are defined as the observations within a contiguous rotation 
range of 180° for the spin axis, and as a consequence each sur¬ 
vey does not cover the full sky. Note that the fifth survey had 
not been completed when HFI stopped observing, and was also 
interrupted by various end-of-mission tests. That last survey was 
also performed with a different scanning strategy than the first 
four surveys (Planck Collaboration I 2015). The date and ring 
number corresponding to the beginning and end of each of the 
time split are summarized in Table A.2. Single survey maps are 
provided for the full channels only. Yearly maps are provided for 
Year 1 and Year 2 for the full channels, the detector sets and the 
SWBs, where the years span Surveys 1-2 and 3-4. Half-mission 
maps are also provided, where each half contains one half of the 
valid rings (or stable pointing periods). There are 347 discarded 
rings (and 26419 valid ones), most of which occurred during 
the (partial) last survey, when various end-of-life tests were per¬ 
formed. 

Half-ring maps are produced by splitting each ring into two 
equal duration parts. The difference of the two half-ring maps 
provides a useful estimate of the high frequency noise and pos¬ 
sibly other systematics. Note that this is the only case where 
the destriping offsets are different from the offsets in the stan¬ 
dard case. Half-ring maps are provided for the full mission only. 
They are given for the full channels, the detector sets, and the 
SWB maps. 

Maps are corrected for zodiacal light emission. Correction 
maps are also provided for the frequency maps and the various 
time splits. 

Units are Kcmb for frequencies up to 353 GHz and 
MJy sr“^ assuming a constant v/y law above. 

A.2. Stokes parameter correlations 

The Planck HFI delivery includes pixelized maps of Stokes co- 
variances {IIJQJU,QQ,QU,UU) solved during the mapmaking 
process for each pixel independently. 

We present in Fig. A.l the distribution of the /, Q, and U 
correlations in each pixel for the HFI frequencies where polar¬ 
ization is reconstructed. 


A.3. Leakage correction maps 

Section 7.3 describes the origin and the formalism of the 
intensity-to-polarization leakage in HFI and how it is, in prin¬ 
ciple, possible to quantify and correct for these systematics ef¬ 
fects. This is particularly important for any studies of the large 
angular scales in polarization, where dipole, monopole, and dust 
leakages are the main limiting factors. 

Although dust and CO bandpass leakage effects can be esti¬ 
mated from given intensity templates and electromagnetic spec¬ 
tra (Sect. A.3.1), such is not the case for the calibration and 
monopole leakage levels, which have to be estimated directly 
from the maps (Sect. A.3.2). 

A.3.1. Bandpass leakage correction maps 

The small differences in the bandpasses of the different bolome¬ 
ters combined to produce polarization maps give rise to some 
leakage from intensity to polarization for the CO and the 
(Galactic) dust emission; this is estimated using the known band 
profiles and templates of the emission (see Sect. 7.3.1). 

Bandpass leakage Q and U maps are provided for dust for 
all channels, using ground-based measurement of the bandpass 
integrated over the dust emission law. Each FITS file contains 
a single extension, with two columns containing the Q and U 
leakage maps in Kcmb units. 

A.3.2. The generalized global fit (GGF) approach for leakage 
correction 

The generalized global fit (GGF) method is a template fitting 
approach that has been developed to consistently solve for cal¬ 
ibration, monopole and bandpass mismatch (BPM) leakage ef¬ 
fects, at the map level. Each Q and U map at 100, 143, 217, and 
353 GHz is modelled as 

[2, ur = [2, U]cmb + [2, (A.l) 

bev bev 

+ F X +F X , 

bev bev 

where: 

- the first line corresponds to the physical polarization signal 
coming from the CMB and the dust (synchrotron is assumed 
negligible at these frequencies); 

- the second line corresponds to the monopole and calibration 
leakage terms, the templates of which make use of the leak¬ 
age pattern maps T described in Sect. 7.3; 

- the dust and CO BPM-induced leakage terms are gathered in 
the third line, again using the leakage patterns maps T; 

- each summation is performed over the polarized bolometers 
of the frequency channel v. 

We use the Planck 353 GHz Q and U maps as polarized dust 
templates and perform the fitting procedure at 1° resolution and 
A^side = 64. Focussing on polarized dust and leakage effects only. 
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Table A.l. Detector set definitions 


Frequency 

DetSetl 

DetSet2 

100 GHz. 

. 100-la/b 100-4a/b 

100-2a/b 100-3a/b 

143 GHz. 

. 143-la/b 143-3a/b 

143-2a/b 143-4a/b 

217 GHz. 

. 217-5a/b 217-7a/b 

217-6a/b 217-8a/b 

353 GHz. 

. 353-3a/b 353-5a/b 

353-4a/b 353-6a/b 


Table A.2. Date and ring numbers for the beginning and end of each time split. 


Time split 

start date 

first ring 

end date 

end ring 

Full mission . 

. 12/08/2009 

240 

13/01/2012 

27005 

Nominal mission. 

. 12/08/2009 

240 

28/11/2010 

14723 

Half mission 1 . 

. 12/08/2009 

240 

15/10/2010 

13471 

Half mission 2 . 

. 15/10/2010 

13472 

13/01/2012 

27005 

Year 1 . 

. 12/08/2009 

240 

12/08/2010 

11194 

Year 2 . 

. 12/08/2010 

11195 

29/07/2011 

21720 

Survey 1. 

. 12/08/2009 

240 

08/02/2010 

5720 
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Eq. (A.l) becomes 
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where is the overall factor scaling the dust from 353 GHz to 
the frequency channel v and where a, JS, y, and S define the am¬ 
plitude of the monopole, dipole, dust and CO leakage effects 
respectively. The last two lines of Eq. (A.2) correct for the to¬ 
tal leakage added when using [Q, as dust templates. This 
equation is at the core of the GGF method, which is then imple¬ 
mented as follows: 

1. A first fit is performed to solve for the coefficients of 
Eq. (A.2), for each channel (v = 100, 143, and 217 GHz) 
independently. There are strong degeneracies between the 
leakage templates at frequency v and at 353 GHz, so that 
this first fit does not provide reliable a, p, y, and S coeffi¬ 
cients that can be used to compute leakage correction maps. 
However, it enables an accurate determination of the over¬ 
all scaling factor of the dust 6^ between the channel under 
scrutiny and 353 GHz. 

2. These values are used as inputs to solve Eq. (A.2), si¬ 
multaneously for all channels, in a consistent manner. At 
this stage, we also add some extra constraints (such as min¬ 
imizing detector-set and survey differences) in order to lift 
further the leakage degeneracies between the four channels. 
This generalized global fit allows us to extract the set of a, p, 
y, and S for all HFI polarized channels, including 353 GHz. 


Solving for survey differences in the second step is crucial. 
Although the leakage coefficients are independent of scan an¬ 
gle, the r maps change because in different surveys pixels are 
scanned at different angles. Total leakage maps, allowing a cor¬ 
rection for monopole, calibration, and BPM leakage, are then 
computed using the coefficients extracted in the second step of 
the procedure. Figure A.2 shows the total correction in Q (top) 
and U (bottom) for the 353 GHz channel. Dust and CO BPM 
leakage effects are clearly dominant near the Galactic plane, 
while the large patterns at high latitude are mainly due to the 
dipole leakage. 
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Fig. A.l. Maps of the correlation IQ (left), lU (middle), QU (right) between Stokes parameters for the 100, 143, 217 and 353 GHz 
(from top to bottom). 
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Fig.A.2. Total leakage correction maps in Q (top) and U (bot¬ 
tom) at 353 GHz computed with the generalized global fit (GGF) 
method. These include calibration, CO and dust leakage tem¬ 
plates. 
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